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' Abstract 

o 

O The classical approach of solvability using group theory is wcill known and one original 

motivation is to solve polynomials by radicals. Radicals are square, cube, square root, cube 

^ root etc of the original coefficients for the polynomial. A polynomial is solvable by radicals 

if the permutation group is solvable. This is exac;t solvability via group theory. With modern 
computers, we might need to relax our definition of exact solvability and move towards practical 

(— I solvability. We will address seismic imaging as an example of practical solvability by semigroup 

^ I theory. The difference between semigroup and group is that the semigroup operators do not 

have to be invertible as in group operators. Using the metaphor of complex decomposition, we 
I will decompose an operator into simple part and complex part. The simple part of the operator 

^ is solvable by numerical methods. The complex part of the operator is interpretable but not 

^sD numerically solvable. It is sometimes called the evanescent energy in geophysics. 

in 

^ Introduction 

Many numerical problems in differential equations and probability can be cast as operators acting 
on Hilbert space. In this paper, we can treat this more abstract operator theory by using n x n 
matrices as a surrogate for the more general operator acting on Hilbert space. Instead of factoring 
polynomials in exact solvability, we can use the same concept of group theory to talk about solving 
the matrix by diagonalization of the matrix with an orthogonal group of matrices. This can be 
viewed as solvability by group theory similar to solving polynomials with permutation groups. If 
we use Gaussian elimination or similar techniques to diagonalize a matrix, each operation on the 
original matrix could be thought of as an invertible transformation of the original matrix. For large 
seismic surveys which are in many terabytes of input data, it is untenable to use the full matrix 
and diagonalize the matrix. One way to generalize the idea is to use semigroup theory where the 
operators are not invertible like rotations in group theory. We can think of taking an operator A 
and write it as sum of idempotent operators, i.e., ^ is a formal sum of PI + P2 + P3 + . . . , where 
A is the operator and Pj is an idempotent with PjPj = Pj- Idempotents are operators P such 
that PP = P. The more familiar notion of idempotent is eigenfunctions or projections or diagonal 
blocks or fixed points or spectral method in numerical analysis. Decomposition into idempotents is 
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a way to achieve stability. We will use wave equation in seismic imaging to demonstrate solvability 
by semigroup theory. 

Seismic imaging 

Migration of seismic data can be classified as operators with predetermined basis and without 
predetermined basis. Migration is the geophysical terminology to map the input recorded data into 
a seismic imaged cube using wave equation. 

The first category includes fk migration and phase shift migration which uses Fourier basis. 
Its advantages are speed and no limitation of bandwidth or dips. The disadvantages are inability 
to handle large velocity variations and inaccuracy of interpolation over a large seismic survey. 

The second category includes finite difference method which includes downward continuation 
and reverse time migration (RTM) and integral method (Kirchhoff or Green's function) . The advan- 
tage is the ability to handle large velocity variations but could have frequency and dip limitations 
due to the practicality of using finite difference over a large survey. Even though technically finite 
difference or Kirchhoff method does not use a predetermined basis, it requires a predetermined 
order of approximation for the whole survey. 

To this end, we will employ an idempotent method which honors large velocity variation and 
all frequencies and dips. It can also accommodate both 1-way and 2-way wave propagation. In 
this paper, we will address the theory and practice of idempotent method. We will not discuss the 
numerical implementation of the method. 

Methodogy 

We will first introduce some mathematical concepts in semigroup and group theory. In order to 
make it less abstract, we will think of operators as matrices since matrices are a good surrogate 
to discuss abstract ideas. If we take all n x n matrices under matrix multiplication as a whole 
system, it is a semigroup and it has the semigroup property of associative law A(BC) = {AB)C. 
This system includes both invertible and noninvertible matrices. We will come back to the idea of 
semigroup and idempotents in semigroups. 

If we restrict a subset to all invertible matrices, then it is a group since each n x n invertible 
matrix A by definition will have a i? so that AB = BA = I the identity matrix. Two matrices are 
equivalent if ^4 = UBV where U and V are invertible. A desirable property is to write a matrix as 
A = UDV where D is a diagonal matrix. More generally, it is desirable to have A = UPV where 
P is an idempotent with PP = P. 

In a group, the only idempotent is the identity. In other words, when we square or cube a 
matrix, it gives a different matrix compared to A. It is deemed desirable to iterate to get a new 
matrix which is different. However, it is in general difficult to know if the iterated matrices converge 
to a meaningful answer. It could converge numerically but it does not mean that it converges to a 
geometrically /geologically meaningful solution. 

An idempotent in a semigroup has been viewed as useless since AA = A. It means that nothing 
new happens. But it is the most stable matrix since nothing new could occur. This is the point 
of using idempotent which is similar to projection operators in numerical implementations. An 
idempotent projects the whole vector space of n-dimensional Euclidean space into a subspace. 

The intuitive notion of an idempotent in migration operator is to select the appropriate pro- 
jection to honor both the wave equation and the complicated velocity model. This is like lumping 
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wave equation and velocity model into ONE operator. The idempotent maintains stability and 
adapts to the large velocity variation (Figure 1). 

A comment on physics versus mathematics 

A typical seismic imaging experiment in oil exploration is to have comprcssional sources (dynamite, 
airguns, vibrating sources etc) at the surface of the earth and gcophoncs/hydrophones to measure 
the energy coming back from the reflections of the earth layers. Our goal is to create a seismic 
cube which can be interpreted as geologic layers to find oil and gas. One such method is downward 
continuation of the wavefield. 

We will use a wave equation without any loss of energy. So for the physics, we will assume 
total conservation of energy. In order to image the earth layers, we need to downward continue 
the wavefield recorded at the surface to deeper layers where the reservoir containing oil and gas 
resides. Downward continuation is the numerical method to extrapolate the wavefield recorded at 
the surface to the deeper layers in the earth. An acoustic or elastic wave equation is used to define 
the extrapolator into the earth. 

Even though the physics is not lossy, the mathematics is lossy. It requires that evanescent 
energy has to be removed at every depth step as we extrapolate the wavefield deeper into the earth. 
Evanescent energy is the numerically unstable part of the recording which grows exponentially in 
each depth step. So even if the physics is conservative, the mathematics has to eliminate part 
of the recorded energy and hence the extrapolation is lossy. It is not a reversible process. If 
we extrapolate into the earth and then undo the extrapolation back to the surface, the wavefield 
without the evanescent energy is not the same as the recorded wavefield. 

The simple part of the wavefield is preserved but the complex part of the wavefield is eliminated 
by the extrapolation method. 

Complex decomposition of operators 

We have used the metaphor of complex decompostion of data into simple part and complex part 
(see Lau et al SEG 2008) . The simple part of the data is numerically explainable and the complex 
part of the data is the left over part which is interpretable but is not explainable by numerical 

methods. This idea was used for data decomposition. With operator decomposition, the simple 
part is the part of the operator without evanescent energy. The complex part is the evanescent 
energy. 

Decomposition of operators (with predetermined basis) 

The simplest way to decompose an operator is to use a predetermined basis like Fourier basis. The 

advantage of predetermined basis is to convert a differential equation into an algebraic equation 
with orthogonal basis like fk (/ is frequency and k is wavenumber). The evanescent wave is just 
filtering certain dips (i.e. certain fk numbers) (Figure 2). The disadvantage is the diflficulty of 
handling large velocity variations. 

Any predetermined basis forces an operator decomposition into predetermined geometry like 
Fourier basis. The predetermined basis maps the operators into fixed shapes given by the basis 
elements. 
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Decomposition of operators (without predetermined basis) 



(1) 



The second method is to lump the derivatives with the velocity field so it is more sensitive to the 
velocity variations. Finite difference and Kirchhoff extrapolation are such methods but they can 
have practical limitations for frequencies and dips. Integral methods like Kirchhoff extrapolation 
(see Berryhill) are more geometrically motivated. An excellent numerical treatment can be found 
in Sandberg and Beylkin 2009 where positive eigenvalues are viewed as evanescent energy. The 
negative and zero eigenvalues are kept as propagating waves in the extrapolation operator. This is 
another way to separate the wavcficld into simple part and complex part. 

We use wave field extrapolation of ID acoustic wave equation to demonstrate decomposition of 
operators without pre-determined basis. ID acoustic wave equation can be written as follows 

^ dz \pdz^ 

where p{z, t) represents pressure, p{z) the bulk density, and c{z) the velocity. A dot above a variable 
denotes differentiation with respect to time. Fourier transform Eq.(l) with respect to time, and 
write in the form of two coupled first-order equations (see e.g., Richards 1971) 



(2) 



where p{z, oj) is the Fourier transform of p{z, t). Equation (2) is written consistently with the con- 
tinuity conditions of continuum mechanics which require that both the pressure and displacement 
remain continuous across all possible interfaces in the medium. 
We can re- write Eq.(2) into matrix form, 
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f = Af, 



which has the following solution when A is independent of z, 

f = e^^fo, 

where io = f{z = 0). 

A traditional way of analysing Eq.(4) is by Taylor series expansion 



(3) 
(4) 



i{z) = ( f; fo = (I + zA + y^A^ + ....)fo 

\fc=o ■ / 



(5) 



Assuming convergence, one can then truncate the series and apply numerical analysis tools to solve 
the equation. 

Another way of further analyzing Eq.(4) is to examine eigenvalues and eigenvectors of A, i.e., 
A = VDV~^, where V is the eigenvector matrix, D diagonal matrix of the eigenvalues. 
Eq.(4) becomes 

f = Ve^"V-^fo, (6) 

where the diagonal eigenvalue matrix is 



and the eigenvector matrix and its inverse are respectively, 
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We can then define a projection of f into g by 

g = V-if. 

Notice the projection depends on the velocity c and the density p. 
Eq.(6) then becomes 

g = e^\o. 

The operator, 

L = e^", 

in Eq.(ll) can be written in more general form as 

L = Pi+P2, 

where 
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(12) 

(13) 
(14) 

(15) 
(16) 

1_ 

are idempotents, which corresponds to upward-propagrating and downward-propagating compo- 
nents of the wavefield, respectively. This simple approach of decomposing ID wave equation oper- 
ator can potentially help in reducing numerical instabilities in wave field downward continuation 
and seismic imaging. 

Topological complexity 

There are geologic layers with fractures and faults which are beyond wave equation or any differ- 
ential equation or probability. This is a more complicated topic since no equation could be written 
when the topology like fractures is too complex. We might need to consider qualitative mathemat- 
ics like computational topology which measures homology or cohomology of the geology. One such 
approach is to compute Betti numbers of the homology or cohomology group. See Kaczynski et al, 
Lau and Yin in arXiv 2010 on L0+L1+L2 optimization. 



5 



Conclusion 



Solvability by semigroup is a way to use noninvertible operators to approximate the original oper- 
ator. We have chosen in this paper to use idempotents P where PP = P. Semigroup theory is not 
restricted to idempotents but idempotents are easier to demonstrate since they resemble diagonal 
matrices (see Davies, Wilansky). They are also familiar to physicists in terms of density matrix 
with or 1 eigenvalues (see Bowler et al) viewed as basic building blocks. 

Multiplication of noninvertible operators becomes less energetic or contractive or lossy or more 
limited. Addition of noninvertible operators becomes more energetic or expansive or generating 
larger operator. Toggling between multiplication and addition of the idempotents gives us the 
stable part (multiplicative) and the generative part (addition). This is the essence of downward 
continuation in wavefield extrapolation. The schematic is captured with the diagram in the end 
(Figure 1). 

The simple part of the wavefield extrapolation operator is written as the sum of idempotent 
operators ( ^ = PI + P2 + P3 + . . . ) where A is the operator and Pj is an idempotent with 
PjPj = Pj and PjPk = if Pj and Pk are different. We dropped the phrase equivalent to idempotent 
to emphasize the importance of idempotents. The simple part of the operator is the propagating 
operator to extrapolate the wavefield. It can be studied through semigroup of operators. The 
complex part of the operator is the evanescent energy which should not be propagated. The 
complex part needs qualitative mathematics like topology to classify and understand. 
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Multiplication Zero Addition 

(less energy) (more energy) 

Figure 1: A lower energy, stable state can be reached by multiplication of idempotents. 
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Input seismic 
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Figure 2: An example of decomposing a 2D seismic section into simple and complex parts, using 
Fourier basis. 



8 



